Phase separation due to quantum mechanical correlations 
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Can phase separation be induced by strong electron correlations? We present a theorem that 
affirmatively answers this question in the Falicov-Kimball model away from half-filling, for any 
dimension. In the ground state the itinerant electrons are spatially separated from the classical 
particles. 



The Falicov-Kimball (FK) model £| can be viewed as 
a modification of the Hubbard (H) model 0] in which one 
species of electrons (say spin down) has infinite mass. As 
such, its relation to the latter is similar to the relation 
of the Ising model to the quantum Heisenberg model || . 
Alternatively, it can be viewed as a model of itinerant 
electrons and immobile ions. It possesses long range or- 
der at low temperature in two or more dimensions at half 
filling, and this checkerboard state (and higher-period 
generalizations PHjJ) remain to-date as the only exam- 
ples of crystallization into a perfectly ordered structure 
whose periodicity is not that of the underlying lattice. 

In this paper, we report a theorem on the existence of 
phase separation in the ground state of the FK model, 
away from half-filling and for large repulsion between 
the particles. 'Phase separation', or 'segregation', means 
that the system splits into two large domains, one being 
occupied by the classical particles, and the other by the 
quantum particles. In the language of the H model, this 
would mean segregation of spin up particles from spin 
down particles — resulting in a ferromagnetic state. The 
question of whether strong interactions can drive quan- 
tum (electronic) systems to phase separate was posed 
over ten years ago for the FK model S and the H model 
H . This work is a rigorous proof of this long-standing 
conjecture for the FK model. Further discussion of the 
relation between the FK and the H model is given later. 

Phase separation is not new in classical lattice mod- 
els. It is present in the Ising model, and in many other 
classical models. It also occurs in the FK model at half- 
filling for some densities, as was proved in l|. In these 
examples it is mainly a local phenomenon — interactions 
(or 'effective interactions' in the case of FK) tend to dis- 
like boundaries, and the state with minimum boundary 
is phase separated. 

Away from half-filling, the electrons of the FK model 
are in delocalized wave functions, and their energy cannot 
be written as a sum of local terms. While we prove that 
the ground state energy is roughly proportional to the 
boundary between occupied and empty sites, the mech- 
anism is nonlocal and genuinely quantum mechanical. 

The FK Hamiltonian iffll is 



x,yGf2 



(1) 



Here, t(x — y) is the hopping coefficient between sites 
x and y; it is translation invariant, but may depend on 
the direction (this allows consideration of general Bra- 
vais lattices). c£ and c x are creation and annihilation 
operators for a spinless electron at site x, and w x = 1 or 
is a classical variable that denotes the presence or the 
absence of an ion at x. (Spin degrees of freedom have 
trivial behavior and are left aside here.) il C Z d is a 
finite d-dimensional lattice, and U is the on-site repul- 
sion between the two species of particles. For any given 
configuration w — {w x } of classical particles, the ground 
state for N e electrons is determined by diagonalizing a 
one-body operator given by the above Hamiltonian, and 
filling in the lowest N c states. The main question is to 
find which configuration w, with a given number of clas- 
sical particles N c = ^ x w x , minimizes the energy of the 
electrons. 

Our theorem states upper and lower bounds for the 
energy of N c electrons, for a given configuration w of 
the classical particles. For orientation, let us consider 
first a configuration where the sites devoid of classical 
particles form a large, 'compact' region. The expected 
energy of the electrons is a bulk term that scales like the 
volume of this region, and a correction that scales like its 
boundary. Our main result is a proof of this conjecture 
for all configurations, not only 'nice' ones. 

We need some notation. Let A = {x G SI : u> x = 0} 
denote the set of empty sites for the configuration w, 
and dA its boundary, <9A = {x G A, dist(x, fi \ A) = 1}. 
Their respective number of sites are |A| and |<9A|. We 
write E(N c ,w) for the ground state energy of N e elec- 
trons in the configuration w. An important quantity is 
n = n e /(l — n c ), where n e = N e /\il\ and n c — ^ x u> x /|S!| 
are the densities for quantum and classical particles re- 
spectively. It represents the electronic density that would 
exist inside A, if all electrons live inside the domain de- 
void of classical particles. Let e(n) be the usual kinetic 
energy per site for noninteracting electrons with density 
n in the thermodynamic limit (its expression is recalled 
below, see Eq. (§)). 
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Theorem. For all A we have upper and lower bounds, 
e(n)\A\ + a'{n)\dA\ > E{N e ,w) > e(n)|A| + a(n, U)\dA\. 

Here, a'(n) and a(n,U) are explicitly given positive 
functions. For nearest-neighbor hoppings (i.e., t(x) ^ 
if |x| = 1, t(x) = otherwise), a(n,U) = a(n) — j(U), 
where a(n) — a(l — n) is strictly positive for < n < 1, 
and 7(Z7) satisfies lim^^oo Uj(U) = 8d 2 . 

The theorem clearly leads to segregation because low- 
energy configurations must have a small boundary, i.e. 
there is a relatively small number of empty sites that 
are neighbors of the classical particles. Furthermore, due 
to the large repulsion, electrons are essentially located 
in the empty domain. The boundary of the configura- 
tion of classical particles that minimizes the energy is 
smaller than the minimum possible boundary times the 
ratio between the upper and lower bounds (which is a 
large constant). 

Phase separation in the FK model for large U was con- 
jectured in |8j and is in stark contrast to the situation at 
n c = n c = 2, where long-range order of checkerboard 
type occurs, as was established in ||. The theorem is 
proved in pj], and extends results for d = 1 and 
d = oo |l2j ] to all dimensions, in particular to the dimen- 
sions 2 and 3 that are of great physical relevance. Its 
proof relies on a result of Li and Yau for the Laplace 
operator in the continuum [jl3| (see also Theorem 
12.3). We explain this elegant proof below in the case 
of a lattice, thereby proving the lower bound for U = oo 
without the boundary correction. To include this bound- 
ary correction requires considerably more effort, and we 
refer to Jl(| for further details. 

Proof with a(n) = and U — oo: The ground-state 
energy e(n) for noninteracting electrons with density n 
is found in the usual fashion: (i) define a Fermi energy 
£f via 



(2tt)° 



d d k = 



e(k)<e F 



(2) 



with the band structure e(k) = — ^ x t(x) exp(— ik • x) 
and (ii) do the integration 



(2ir) d 



£(k)<e F 



e(k) d d k = e{n). 



(3) 



The electrons are forbidden to lie on any site occupied 
by the classical particles. Then the eigenfunctions for 
a given configuration w are found by diagonalizing the 
projection of the hopping matrix onto A. Let 4>f3(y;w) 
denote the orthonormal eigenvectors of the Hamiltonian 
indexed by (3 — 1, |0| — N c for the configuration w 
with eigenvalues efl(w). We choose the ordering of the 
labels such that ei(w) < £2(10) < ••■ < £|n|-JV c i w ). We 
also extend the definition of the eigenvectors to all of Z d 
setting <t>p(y;w) — for all y ^ fl. The ground state 
energy for iV c electrons and iV c classical particles in the 
configuration w is 



13=1 



(4) 



and the ground-state energy E g , s ,(N c , N c ) is the mini- 
mum of Eq. (jfy over all configurations that contain 
classical particles. 

Using the definition of the eigenvectors allows us to 
write the ground state energy as 

£(tv , w ) = ^ J2 ^(y; w )H(y- z )]^( z ;H- (5) 

13=1 y,z£Z d 

Now, we define the Fourier transform of the eigenfunc- 
tions by 



Mk;w)= £ e ik ^( yri 



(6) 



for every wave vector k in the Brillouin zone [0, 27r] d . It 
is well-known that the energy can be expressed as 



E(N c ,w) = ^ J d d ke(k)p(k;w); 



(7) 



here and in the following the integral is over [0, 27r] d and 
we introduced the density function 



p(k; W )=^|/ /3 (k;«;)| 5 



The density function is obviously positive and is bounded 
above by |A|. Indeed, we can write 



p(k;w) = Y. E e-^&wy^Mv-,™)-, (9) 

0=1 y,z6Z d 



this can be rewritten as 



p(k;w)= ]T e- ik >(z, y;U ;)e ik y 

y,z£A 



(10) 



with p(z,y;w) = ^9=1 ^( z ! w )^(y'i w )' The matrix p 
is a positive semidefinite matrix bounded by 1 (it is the 
projector onto the lowest A^ eigenvectors). Hence, 



p(k;u,)<5> ik ->f = |A| 

yeA 



(11) 



Furthermore, the density function satisfies a sum rule 

-L I d d k P (k: w) = J2 p(y» y; w) = N e . (12) 

"" 1 " yeA 



(27 



One gets a lower bound for E(N C , w) by minimizing the 
rig ht side of (|) over all functions p satisfying < p < \A\ 
and whose integral is AT C . This is the 'bathtub principle', 
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see [jlij, Theorem 1.14; the minimizer is p(k) = |A| for 
{k : e(k) < £f} and p(k) = otherwise, with ep defined 
by Eq. (||). This amounts to filling the lowest eigenvalues 
of the infinite-lattice. 

The proof of the upper bound for the energy proceeds 
by forming an 'average' Hamiltonian by translating and 
rotating the original configuration w over a large but fi- 
nite subset of Z d (with periodic boundary conditions). 
Then on a bipartite lattice, one can show by concavity 
of the sum of the lowest N eigenvalues of a matrix, that 
the averaged Hamiltonian provides an upper bound to 
the ground state energy. But the magnitude of the aver- 
aged hopping is determined by the size of the boundary, 
which eventually yields the desired upper bound. Ex- 
tending the proof above to provide the lower bound is 
much more complicated and relies on a detailed techni- 
cal examination of the influence of the boundary sites on 
the minimal density function for a given configuration of 
classical particles. It is done in the isotropic case (that 
is, a hypercubic Bravais lattice with equal hoppings in all 
directions) in |l(| ; the extension to the anisotropic case is 
straightforward, but numerous details are modified, and 
have been verified. 

The results of the theorem have a number of impli- 
cations for the FK model. It establishes that the seg- 
regation principle holds in all dimensions (at T = and 
U = oo) illustrating the fact that the existence of periodic 
ground states requires a subtle reduction in energy rela- 
tive to the segregated phase as the interaction strength 
is made finite. Since the electronic wave functions will be 
exponentially localized within A, the results shown here 
can be extended to the case of finite interaction strength, 
as long as U is large enough. At positive temperature it 
is so far impossible to claim rigorous results, except the 
following weaker one: the electronic free energy (for a 
fixed configuration of classical particles) can be shown to 
be equal to the bulk free energy plus a correction term 
that is proportional to the size of the boundary. We ex- 
pect that the coexistence of two phases occurs at finite 
temperature for d > 2 as happens in the Ising model. 

It may be instructive to consider the U — oo results 
for the FK model as a guide for possible behavior in the 
H model. To do this, we must first find a way to in- 
terpolate between the two models. The simplest way is 
to consider an asymmetric-hopping H model where the 
hopping for the spin-up and the spin-down particles is 
different. Then the H model results when if = if and 
the FK model when if = 0. Our rigorous results only 
hold for ty = 0. When if is increased the classical par- 
ticles should still be packed, due to the pressure of the 
electrons. When if keeps increasing however, the clas- 
sical particles should be in a phase with density strictly 
less than 1. The central issue is whether the reduction of 
the down-spin kinetic energy can be made large enough, 
so that the phase separation disappears at a critical value 
of if . If this occurs for all electron densities, there is no 



(saturated) ferromagnetism in the U = oo H model on 
the given lattice; however, one has a saturated ferromag- 
netic ground state if the phase separation survives (note 
the SU (2)-imposed degeneracy of the ferromagnetic mul- 
tiplet will occur precisely at if = if). It is well known 
that ferromagnetism depends strongly on the geometry 
of the lattice |l5|,[l6| , so the occurrence of a critical value 
of if must also depend strongly on the geometry of the 
lattice. We are unable to make any rigorous statements 
about ferromagnetism in the H model here — actually, 
we do not even know how to study the case with nonzero, 
but small if. 

A major question is what happens to the chessboard 
phase when doped away from half filling? Consider the 
line n e = n c ; the chessboard phase is present when these 
densities are equal to \ , and segregation takes place when 
they differ significantly from ^ (depending on U). It is 
not clear what to expect for intermediate values. Two 
possible scenarios are (i) the coexistence between chess- 
board and segregated phases or (ii) the coexistence be- 
tween other periodic and segregated phases. Both sce- 
narios could be of physical relevance to stripe physics. 

We conclude this letter by a summary of our knowledge 
of the phase diagram for zero temperature and d < 2. 
Recall that the particle-hole symmetry implies that the 
phase diagram is symmetric under the transformation 
(n c , n c ) (1 — n c , 1 — n c ) M. In the sequel we describe 
the situation for n c + n c < 1, which is enough. 

In two dimensions, the ground states are periodic 

• when (n e , n c ) = (i, |) (they are of the chessboard 
type). This was proved for all U in ||; 

• when n c = 1 — n c (i.e. at half-filling) , and n c — 



2 1 1 2 1 _2_ 1 

5' 3' 4' 9' 5' 11' 6 



also, n c = 



with intc- 



n 2 + (n+l) 2 

gers n. This holds for U large enough (depending 
on n c ), and follows from |g-|,||,Q; 

• when n c = 1 — n e , and the electronic density n e is a 
rational number between | and | . U must be larger 
than a value that depends on the denominator of 



There is coexistence of two periodic phases 

• when n c = l-n e and n e € (I, JL) y (|, §) U (§, ±), 
for U large 



And the ground states display segregation 



• for 



Hi 



< n c < (1 



described here, and proved in 



)(1 — n c ); this is 



These results are illustrated in Fig. [l| 

The domain c ™ st (1 — n c ) > n c should also be segre- 
gated. The central band that includes the line n c +n c = 1 
should be the host of numerous periodic phases and var- 
ious coexistences between periodic phases and empty or 
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at Princeton University where the majority of this work 
was completed. We thank S. Kivelson, R. Lemahski and 
M. Loss for helpful discussions. 



density of electrons 

FIG. 1. Schematic phase diagram of the rigorous results for 
the ground state of the 2D Falicov-Kimball model for large 
U. The white dots and the white lines represent periodic 
phases; the black lines are coexistences between different pe- 
riodic phases; and dark gray regions are segregated. There 
are no rigorous results for light gray domains. 



full phases. This is supported by numerical simulations 
in 2D @. 

In one dimension, the ground states are periodic 



when (n c ,n c ) = (|, ±) 



• when n c = 1 — n c (half-filling) and n c — p/q is a 
rational number (in an irreducible fraction). The 
periodicity is q for U sufficiently large (depending 
on?) § 

For finite U, there is numerical evidence for coexis- 
tence of two periodic phases and free electrons 

• when n c = 1 — n e , p/q < n c < p' /q', and the peri- 
odic phases with period q and q' are the only stable 
phases within the above interval pfj[ . 

And the ground states display segregation 

• when n e 1 — n c and U is sufficiently large pTT^pToj . 

The canonical phase diagram for small U is even richer 
but our knowledge of it is very limited. 

In conclusion, we have proved that the FK model is 
phase separated whenever N c < |A| and U — > oo. This 
shows how strong correlations can lead to phase separa- 
tion. 
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